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Using the method of maximal cuts, we construct the complete D-dimensional integrand of the 
five- loop four-point amplitude of Af — 4 super- Yang-Mills theory, including nonplanar contributions. 
In the critical dimension where this amplitude becomes ultraviolet divergent, we present a compact 
explicit expression for the nonvanishing ultraviolet divergence in terms of three vacuum integrals. 
This construction provides a crucial step towards obtaining the corresponding amplitude of M — 8 
supergravity required to resolve the general ultraviolet behavior of supergravity theories. 
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PACS numbers: 04.65.+e, ll.15.Bt, 11.30.Pb, 11.55.Bq 

Recent years have seen remarkable progress in under- 
standing and constructing scattering amplitudes in gauge 
and gravity theories, driven largely by the advent of on- 
shell techniques. The advances have had broad applica- 
tions including QCD computations of multijet processes 
at the Large Hadron Collider, resummations of planar 
Af = 4 super- Yang-Mills (sYM) amplitudes linking them 
to string theory via the AdS/CFT correspondence, con- 
nections to integrability of planar J\f = 4 sYM theory and 
studies of ultraviolet (UV) and infrared diverges in gauge 
and gravity theories. (See ref. [l[ for recent reviews.) 

These advances have been most striking for the max- 
imally supersymmetric Af = 4 sYM amplitude, in the 
planar limit where the number of color charges is large. 
Significant progress has also been made for the less well 
understood nonplanar case, which is the subject of this 
Letter. Nonplanar contributions to amplitudes in Af = 4 
sYM theory have been obtained previously through four 
loops [2-7] , along with detailed studies of their UV prop- 
erties in higher space-time dimensions. The planar part 
of the five-loop amplitude is found in ref. [8j. Here we 
carry out a similar study for the five-loop four-point am- 
plitude and analyze the UV divergences in the lowest di- 
mension where they occur. We anticipate that our results 
will become useful for detailed studies of the structure 
of the theory, including infrared singularities, anomalous 
dimensions and other observables related to amplitudes. 

Beyond the intrinsic interest for understanding Af = 4 
sYM theory, our construction of the five-loop four-point 
amplitude is a key step towards obtaining the correspond- 
ing amplitudes of Af > 4 supergravity, needed to help 
resolve the long-standing question on the possible UV 
finiteness these theories. In fact, whenever a represen- 
tation of an Af — 4 sYM amplitude is constructed that 
exhibits a duality between color and kinematics 
simple pathway exists for obtaining corresponding Af > 4 



supergravity loop integrands @, 0, El, El- In particu- 
lar, the Af — 8 supergravity integrand follows trivially 
when the duality is manifest, simply by replacing color 
factors with the kinematic numerators of the diagrams. 
Although the form of the five-loop four-point amplitude 
presented here does not manifest the required duality, 
it does offer an excellent starting point for finding such 
representations. 

Explicit constructions of amplitudes have played a key 
role for determining the UV divergence structure of gauge 
and gravity theories as a function of dimension. Af = 4 
sYM theory is known 0, [l3[ to be UV finite in dimensions 
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where L is the loop order. This exhibits the well known 
UV finiteness in D = 4 [14]. A remaining open question 
is whether the bound fl} is saturated to all loop orders. 
Prom explicit computations, it is known to be saturated 
for L < 4 [H 0, 0, 0] ■ As commented on in ref. @ , the 
L = 5 planar amplitude [1] also saturates the bound ((T|). 
Below we give a simple expression for the divergence, 
including nonplanar parts. 

A related open question is whether maximally super- 
symmetric Af = 8 supergravity has the same finiteness 
bound as Af — 4 sYM theory, implying it is UV finite 
in D = 4, or if it has a worse behavior. (For a re- 
cent optimistic opinion see ref. [l5| ; for a recent pes- 
simistic one see ref. I 16fl) Explicit calculations of the 
divergences 0, 0, 0, Il7l Il8j and symmetry and other ar- 
guments [3] show that through four loops the bound ([T]) 
holds in Af = 8 supergravity. Beyond this, the argu- 
ments suggest that Af = 8 supergravity will have a worse 
behavior, leading to a seven-loop divergence in D = 4. 
However, when similar symmetry arguments are applied 
to Af = 4 supergravity, they imply the existence of a valid 
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FIG. 1: Sample graphs for the five-loop four-point M = 4 
sYM amplitude. The graph labels correspond to the ones 
used in the ancillary file [231 ]. 



three-loop counterterm [20(; the coefficient of this coun- 
terterm has recently been explicitly shown to vanish [l2| . 
(See ref. 21| for a string-based argument.) This exhibits 
better behavior than implied by known symmetry con- 
siderations and is in line with cancellations suggested by 
unitarity arguments 22]. In particular, it emphasizes the 
importance of directly checking the amplitudes whether 
cq. (Q]) holds for TV = 8 supergravity at L = 5. 

Our construction of the five-loop four-point amplitude 
of AT = 4 sYM theory organizes it in the form, 
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where the second sum runs over a set of 416 distinct 
(non-isomorphic) graphs with only cubic (trivalent) ver- 
tices. Some sample graphs are shown in fig. [1] The 
first sum runs over all 24 permutations of external leg 
labels indicated by £4. The symmetry factors Si remove 
overcounts, including those arising from internal auto- 
morphism symmetries with external legs fixed. Here we 
absorb all contact terms (i.e. terms with fewer than the 
maximum number of propagators) into graphs with only 
cubic vertices, by multiplying and dividing by appropri- 
ate propagators. We denote external momenta by ki for 
i = 1, . . . , 4 and the five independent loop momenta by 



j for j = 5, . 



,9. The remaining lj are linear combi- 



nations of these. The color factors Ci of all graphs are 
obtained by dressing every three- vertex in the graph with 
a factor of f abc = Tr([T a , T b ]T c ), where the gauge group 
generators T a are normalized as Tr(T a T b ) = S ab . The 
gauge coupling is g and the crossing symmetric prefac- 
tor stAf co is in terms of the color-ordered D-dimensional 
tree amplitude A\ rcc = A\ rcc (l, 2, 3, 4) and s = (k 1 + k 2 ) 2 
and t = (fc 2 + fc 3 ) 2 . 

To construct the numerators Ni , we use the method of 
maximal cuts Q , based on the unitarity method 24 1 . Ap- 
plication of this method and various strategies for greatly 
streamlining the construction of the numerators has been 
described in considerable detail in ref. @, so here we 
give only a brief summary. The method works in D di- 
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FIG. 2: Sample N fc -maximal cuts for k = 0,1,2,3. The ex- 
posed lines are all cut. 
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FIG. 3: Examples of simple cuts used to speed up the calcu- 
lation, (a) is a two particle cut, (b) a box cut and (c) is a 
sample application of BCJ amplitude relations. The exposed 
lines are all cut. 



mensions and can be used to obtain local expressions, 
from which UV divergences can be straightforwardly ex- 
tracted. 

We start with an ansatz for the diagram numerators 
containing free parameters to be determined by matching 
against generalized unitarity cuts. Our ansatz is a poly- 
nomial of degree four in the kinematic invariants, subject 
to the power-counting constraint that no term has more 
than six powers of loop momentum. We also demand 
that each numerator respects the automorphism symme- 
tries of the graph. Once a solution is found satisfying 
a complete set of cut conditions, we have the integrand. 
If an inconsistency is encountered, the ansatz must be 
enlarged. We note that the solutions for numerators are 
not unique and different choices can be mapped into each 
other by generalized gauge transformations 

The parameters of the ansatz are determined from gen- 
eralized unitarity cuts that decompose a loop integrand 
into products of on-shell tree amplitudes summed over 



all intermediate states, £ ?tates AfifAfi* ■ • • Ap}. These 
cuts are organized according to the number of cut prop- 
agators that are replaced with on-shell conditions. We 
start from the maximal cuts (MCs) where all 16 internal 
propagators cut. After obtaining the MCs, we then con- 
structs all next-to- maximal cuts (NMCs), with 15 cut 
propagators. We continue this process, systematically 
constructing analytic expressions for (next-to) ^-maximal 
cuts (N fc MCs) with fewer and fewer imposed cut condi- 
tions. For the five-loop four-point M = 4 sYM amplitude 
this process terminates at k = 3, since the power count- 
ing of the theory prevents numerator factors from can- 



celing more than 3 propagators. Representative cuts for 
k = 0, 1, 2, 3 are shown in fig. [2] The number of nonzero 
(color-stripped) cuts of type N fe MC arc 410, 2473, 7917, 
15156 for k = 0,1,2,3, respectively. This count does 
not include the different independent color orderings of 
each cut. In addition to the nonvanishing cuts, there is 
a large class of N fc - 3 MCs that evaluate to zero because 
they contain nontrivial (n < 3)-point subamplitudcs. 

Each cut can be reduced to a relatively simple analytic 
expression. The N fc MCs are evaluated in D dimensions 
by embedding them in auxiliary cuts that can be directly 
expressed in terms of simplified analytic forms. As dis- 
cussed in some detail in ref. [g], two particularly useful 
cuts for this purpose are two-particle cuts and box cuts. 
Whenever a two-particle reducible cut can be factorized 
into two four-point amplitudes, as illustrated in fig-EJa), 
all contributions to the cut can be written down immedi- 
ately using lower-loop results [2j. Similarly, all cut con- 
tributions that possess a four-point loop or box subdia- 
gram, illustrated in fig. G^b), are simple to evaluate [6j. 
A third type of auxiliary cut 9] , illustrated in fig. ^c) , 
allows us to map known D-dimensional planar cuts to 
nonplanar ones via the BCJ tree-amplitude relations 
Alternatively, one can construct numerator representa- 
tions that obey the color-kinematics duality for each cut 
separately [6j, giving local numerator relations between 
planar and nonplanar diagrams, up to terms that vanish 
on the cut. This technique is especially useful whenever 
the cut contains massless bubble or tadpole subdiagrams 
(as sometimes occurs for N 3 MCs), since the local numer- 
ators are automatically free of spurious singularities that 
can appear with other methods. A fourth type of auxil- 
iary cut used in our construction is one that maps five- 
loop nonplanar cuts to simpler six-loop planar cuts j2(| . 

We have found a choice of parameters in the starting 
ansatz whose cuts correctly reproduce the N fc MCs at the 
level of the integrand. We thus have a complete integral 
representation of the five-loop amplitude. As a few sim- 
ple examples, the numerators of graphs 1, 12 and 284 
are 

Ni = s 4 , A^i 2 = 2s 3 fc 3 • l 5 , 

A 284 = 2 S 2 ((/ 10 -/ 2 o) 2 + ai3-/i8) 2 ), 

corresponding to the graphs in fig. Q] labeled as (1), (12) 
and (284) and matching the labeling in the ancillary 
file [231 ] . The lines with arrows in fig. [1] give the momen- 
tum labels and directions. The symmetry factors 1/Si 
for these graphs are respectively 1/4, 1 and 1/4. 

The complete set of 416 nonvanishing graphs with their 
associated symmetry factors, numerators and color fac- 
tors are included in the ancillary file 23]. We note that 
graphs 61, 67, 133, 137, 263, 382, 412 have vanishing color 
factors for a general gauge group (due to symmetry prop- 
erties of the graph), and hence do not contribute to the 
amplitude. However, we include them in our represen- 
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FIG. 4: Some of the five-loop vacuum integrals that appear in 
intermediate steps. Only (a), (b) and (c) appear in the final 
UV divergence. Integral (j) has a nontrivial numerator factor, 
as indicated. The (blue) dots indicate that a propagator is 
squared. 

tation because they are needed for constructing gravity 
amplitudes [6|. 

We have carried out extensive cross checks on our re- 
sult. The cut construction automatically cross checks 
the vast majority of contributions because they are de- 
tected in multiple independent channels. As an addi- 
tional rather nontrivial check, in four dimensions we con- 
firmed numerically that all the analytically-obtained cuts 
are correct; to carry out this check we used the simple 
algorithms of ref. [27j for carrying out the supersums ap- 
pearing in the cuts. We have also carried out systematic 
cross checks using generalized cuts with up to six col- 
lapsed propagators. Furthermore, we have evaluated a 
set of cuts that suffices to detect all "snail" contribu- 
tions, equivalent to bubbles on external legs (see sections 
2D and 3C of ref. Q), showing that such contributions 
do not appear in our representation. 

Starting with the constructed integrand, we obtain the 
potential logarithmic divergence in the five-loop critical 
dimension, D = 26/5, following the same strategy as at 
lower loops 0, 0, E3, [Hj] : We expand the amplitudes 
at small external momentum and keep the leading term. 
The result of this expansion is a sum of about 185 vac- 
uum diagrams; a few of which are displayed in fig. [4] 
As discussed in refs. 0, EH, the vacuum integrals in this 
expansion are not all independent (so the precise num- 
ber appearing initially can vary). We derive consistency 
relations between the vacuum integrals by considering 
auxiliary linearly divergent integrals of similar propaga- 
tor structure, expanding them around zero external mo- 
mentum and requiring that the results of the expansion 
be independent of different integrand parametrizations. 
This also directly cross checks the procedure for inte- 
gral reduction since we obtain a highly overconstrained 
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set of homogeneous consistency equations. The fact that 
no positive definite integral is set to zero by this system 
is a strong check on the calculation. These consistency 
relations eliminate most of the vacuum diagrams. Two 
examples are, 
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where the labels correspond to the ones in fig. 2J 

After using the consistency relations, the leading UV 
divergence is remarkably simple and given by only three 
vacuum integrals. For SU(N C ), it is 
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With the chosen normalization, the Wick rotated vacuum 
integrals in eq. ([3]) are all positive definite, proving that 
no further hidden cancellations remain at L = 5 in the 
critical dimension for either leading- or subleading-color 
contributions. Using FIESTA [28[ we have numerically 
evaluated the integrals giving, 
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where the dimensional regularization parameter is e = 
(26/5 — D)/2, K = l/(47r) 13 and numerical integration 
uncertainties are below the displayed digits. It is inter- 
esting that the ratio between the subleading and leading 
contributions 45.0/7V^ is rather close to the three- and 
four-loop ratios, 43.3/Y c 2 and 44.4/Y c 2 Q. A striking 
feature of the result (O is that the divergence does not 
contain terms beyond 0(1/ N 2 ) suppression, nor does it 
contain double-trace contributions when converted to an 
SU(N C ) color-trace representation, in line with expecta- 
tions from lower loops The second of these features 
has already been discussed in refs. [l,[29j]. Furthermore, 
the three integrals and their relative coefficients have a 
remarkable similarity with the corresponding ones at four 
loops, as seen by comparing to eq. (5.33) of ref. [fj. At 
lower loops, exactly the same combination of integrals 
appearing in the subleading-color contributions to the 
Af = 4 sYM divergences appear in the corresponding 
ones of M = 8 supergravity [6j. A natural conjecture is 
that the same holds at five loops, so that the two theories 
share the same critical dimension, D = 26/5. 

In summary, the five-loop amplitude we have con- 
structed here offers detailed information on the struc- 
ture of the nonplanar sector of J\f = 4 sYM theory. As 



a first application, we have shown that simple patterns 
for divergences in the dimension where they first appear 
continue to hold through five loops; this hints that the di- 
vergences are controlled by a deep structure of the theory. 
Our construction of the five-loop four-point amplitude is 
an excellent starting point to try to find a representation 
exhibiting the duality between color and kinematics. We 
expect that the results presented here will be crucial in- 
put for obtaining corresponding supergravity amplitudes 
and for studying their UV behavior. 
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